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ABSTRACT 

The 7-rays from 7-ray bursts (GRBs) are believed to be produced by internal shocks driven by small 
timescale, ~ 1 ms, variation in the GRB outflows, and a pair-production spectral cutoff is generally expected 
around the GeV range. However, the observed optical flashes accompanying GRBs suggest that the delayed 
residual collisions due to large timescale variation continue to accelerate electrons. We show here that the 
inverse-Compton (IC) scattering of the prompt 7-rays by these residual internal shock electrons leads to a high 
energy emission beyond the previously thought spectral cutoff, in agreement with the previous detections of 
GeV photons by EGRET in several GRBs in conjunction with MeV emission. We expect a spectral break 
due to the transition from the primary to residual internal shock emission at the previously thought spectral 
cutoff, and expect systematic time delays of high energy photons relative to MeV emission, the discovery of 
which would provide stringent constraint on the outflow properties, but requires large enough collection of high 
energy photons by, e.g., Fermi and AGILE satellites. 

The recent Fermi-detected bright GRB 080916c unambiguously shows the shifting of the prompt emission 
toward later times as the photon energy increases. The second-scale shifting at >100 MeV is much longer than 
the MeV variability time, as predicted in the residual collision model. The observations imply that there should 
be emission above 70 GeV in the source frame, which may not be produced by primary internal shocks but by 
IC emission in residual collisions. With the method involving time delays of high energy emission, the bulk 
Lorentz factor of GRB 080916c is determined to be L 300. 

Subject headings: acceleration of particles — magnetic fields — shock waves — gamma-rays: bursts 
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1. EVITRODUCTION 

The prompt MeV 7-rays from a GRB are well believed to 
be produced by a unsteady outflow which causes internal col- 
lisions between different parts with various velocities, leading 
to kinetic energy dissipation (Paczynski & Xu 1994; Rees & 
Meszaros 1994)(see, however, Lyutikov & Blandford 2003; 
Narayan & Kumar 2008). The internal shock model can nat- 
urally explain both the non-thermal spectra and the compli- 
cated light curves of GRBs. The observed temporal variabili- 
ties of 7-ray emission are believed to be reflecting the activi- 
ties of the central engine (Sari & Piran 1997; Kobayashi et al. 
1997). The internal shocks are expected to generate/amplify 
magnetic field and accelerate electrons, leading to MeV 7- 
rays by synchrotron emission (see, e.g. Waxman 2003, for a 
review). 

GRBs are mainly observed in MeV range, the properties 
of high energy, say, > lOOMeV, emission are not well under- 
stood, which, however, might be very helpful in constrain- 
ing the physics of the GRB emission region. For example, 
the observed > 100 MeV photons in conjunction with MeV 
emission by EGRET in several GRBs, suggesting that they 
can avoid the 77 absorption, have leaded to the conclusion 
that the GRB outflow must be relativistically expanding with 
a Lorentz factor of F > 100 (e.g. Baring & Harding 1997)(see 
also Waxman 2003). As the development in high energy 
7-ray observations, e.g. AGILE and Fermi are being well 
operated, there are great interests on detecting the high en- 
ergy pair-production spectral cutoff in order to constrain the 
size and Lorentz factor of GRB emission region (e.g. Baring 
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(2000); Lithwick & Sari (2001); and recent detailed consider- 
ation by Gupta & Zhang (2007, 2008); Granot et al. (2008); 
Murase & loka (2008)). 

However, a simple spectral cutoff may not exist. Li & Wax- 
man (2008, hereafter LW08) had noticed that the frequently 
observed prompt optical emission (Akerlof et al. 1999; Blake 
et al. 2005; Vestrand et al. 2005, 2006; Yost et al. 2007) 
is above the expected synchrotron-self absorption limit, and 
suggests a relatively large size of optical emission region, 
compared to that of MeV emission. Actually in the context 
of internal shock models one would expect delayed collisions 
in the outflow following the small radius collisions driven by 
smallest timescale variation of the outflow properties, and 
these delayed collisions can naturally account for the rela- 
tively bright optical emission (LW08). The recently detected 
optical flash from the naked-eye GRB 080319b (Bloom et al. 
2008; D'Elia et al. 2008; Racusin et al. 2008; Wozniak et al. 
2008) appears to vary rapidly in times and its temporal profile 
is correlated to the MeV emission in second scales (Beskin 
et al. 2009), supporting that the prompt optical emission in 
GRBs is produced by internal shocks within the outflow, i.e. 
synchrotron emission from residual collisions (Li & Waxman, 
in preparation). The electrons in residual collisions mainly 
cool by IC scattering off the MeV photons, which produces 
high energy emission at large radii where the optical depth due 
to pair production is reduced (LW08). We consider here this 
high energy emission and show that it may "smear" the previ- 
ously thought pair-production spectral cutoff. But a spectral 
turnover is still expected, which may be observed by Fermi 
and AGILE although more difficult to detect than a simple 
cutoff. The high energy emission from residual collisions is 
also expected to be delayed relative to MeV 7-rays . It should 
be noticed that we only focus on the prompt high energy emis- 
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sion which appears simultaneously with the MeV 7-ray emis- 
sion. 

We show first in §2 the strong 77 absorption during the 
MeV 7-ray emitting phase, next discuss in §3 the high en- 
ergy emission from residual collisions, then the model is ap- 
plied to the recent Fermi detection of GRB 080916c (Abdo et 
al. 2009) in §4, which might have provided evidences of the 
model, finally a general discussion on observations is given in 
§5. 

2. 77 ABSORPTION AT SMALL RADII 

Consider a highly relativistic outflow with a bulk Lorentz 
factor r = 10^^r2.5. The small timescale variation will lead to 
strong internal shocks, which produce 7-ray emission. Let us 
denote the radius where 7-rays arise by R^. Due to geometry 
effect, the observed fastest variabihty timescale ?var ^ I0~^s 
(e.g. Woods & Loeb 1995) in 7-ray light curves suggests that 
the size of 7-ray emission region is limited to < 2r^cfvar- 
If the 7-ray emitting electrons are fast cooling (with cool- 
ing time shorter than the dynamical time of the outflow), 
we should take the equality, thus, /?-y w 6 x lO^^Fj j^var -2cm, 
where fvar = 10~^fvai-.-2S. Actually, in the context of inter- 
nal shock models, we do not expect the magnetic field is 
generated with energy density much higher than that of ac- 
celerated electrons. In order to have synchrotron emission 
peaking at ~ 1 MeV, as observed in GRBs, the radius 
of the 7-ray emitting region should be small (LW08), R^ < 

10'^L^^j^52(efc/l MeV)"'cm, where Lboi = 10^%oi,52erg s"' is 
the bolometric 7-ray luminosity. 

Let us consider the 77 absorption due to pair production 
inside the GRB source. For a photon of high energy e the op- 
tical depth to pair-production within the GRB source is given 
by the product of the pair-production rate, l/t!y^{e), and the 
dynamical time, the time required for significant expansion of 
the plasma, t'^ ~ R^/Tc (primes denote quantities measured in 
the outflow rest frame), T~^j(e) ~ R~f/Tct!^^{e). f '^,^(e) depends 
on the energy density and on the spectrum of the radiation. 
The (outflow rest frame) radiation energy density is approx- 
imately given by t/^ = L/4'!TR^cr^. The GRB spectrum can 
typically be described as a broken power law, dn /de oc e"'', 
with /3 f» -1 at low energy, e < ~ 1 MeV, and /3 « -2 at 
€> eh (Band et al. 1993). High energy photons with energy 
e' exceeding e'/^ = linieC^)^ / e'l^, may produce pairs in interac- 
tions with photons of energy exceeding e' = KnieC^Y / e' < 
(the rest frame photon energy e' is related to the observed en- 
ergy by e = Fe') • For e' < we have dn/de' (xe'~^ , which im- 
plies that the number density of photons with energy exceed- 
ing e' depends only weakly on energy. Thus, t^^ is nearly in- 
dependent of energy for e' > t'~^ ~ {aj / 16)ct/^ /e^, which 
gives 



T^^(£ > Si) ~ 1.1 X 10" 
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fvar,-2ri5 VlMeVy' 



(1) 



Note, we have approximated the 77 cross section above 
the pair-production threshold as 3(Jr/16 (Waxman 2003). 
Also note that as the energy density U'^ (and hence the re- 
lated luminosity L) depends on the energy band considered, 
the one used in calculating t'^^ is the energy density below 
2 X e'ly Hereafter, without special emphasis the luminosity 
L = 10^^L52ergs~' is the so-called MeV luminosity, only cor- 
responding to emission at < 2 x ej,, i.e., L = J^^'' L^de. 



Photons of lower energy, e < Sh, interact to produce pairs 
only with photons of energy e' > 2(meC^)^/e' > e^. Since 
the number density of these photons drops like 1 / e', we have 
Tj^ie < Sb) w i£/eb)T^jie > £&), i.e. 



T-y^(e < eb) ^ 2.2 x 10 
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The optical depth increases as photon energy increases. Pho- 
tons with high enough energy are absorbed in the emission 
region. A spectral cutoff is defined by r-y^ = 1,^ 



e^, ~ 0.46 



GeV. 



L52 



(3) 



A comment on the approximation of the 77 cross section 
should be made here. Since both the cross section and the 
GRB photon spectrum decrease rapidly with photon energy, 
the approximation is excellent- for a GRB spectrum with 
P « —2, it precisely produces the optical depth within 2%, 
compared to a calculation with full cross section. 

There is another restriction for the cutoff besides eq. (3). 
In deriving the cutoff, eq. (3), the optical depth is not self- 
consistently calculated since a GRB spectrum extending to 
infinity without a high energy cutoff is assumed. Given the 
two factors that GRB spectrum usuaUy appears to be a steep 
slope, with the photon number dominated by low energy pho- 
tons, and that the high (e' > l^l^meC^) and low (e' < l^l^meC^) 
energy photons annihilate each other one by one, we only ex- 
pect high energy photons might be totally attenuated by low 
energy ones, other than the other way around. So the cutoff 
should be e^^ > which is not assured by e^ul in eq. 

(3) (Scut/r < 2^l^meC^ might happen). We need to define, in 
the GRB source frame. 



efl = l^l^Tm.c^ ~ 0.23r2.5 GeV. 



(4) 



Note, this condition is also forgotten by many other authors 
who calculate the cutoff energy assuming a no-cutoff GRB 
spectrum. The GRB spectral cutoff is, therefore, expected at 



the maximum between and e^ut (Li & Song 2004), 



(2) 



£0 = ^catiR-i) — max 



0.46 



fvar,-2F2 5 
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O.23F2.5 



GeV. (5) 



A critical Lorentz factor where e*J,t = ^cut is function of R-^ 



(hence ?var). 



:280(L52Avar.-2)''''. 



(6) 



£0 sensitively depends on F thus detection of the spectral cut- 
off may be very useful in constraining F, 



F = min 
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t°° 1 GeV 

, 'var,-2 i V 



(l + z)i/3, 1390^^(1+2) 
(7)' 

Here the redshift factor has been included, i.e. = t{l +z) 

and £°'^ = £/(l+z). 

The extension of GRB spectra to > 100 MeV and the char- 
acteristic variability time, ?var-2 ^ 1, have implied r2.5 > 1, 

^ Because a low energy turnover at ~ 1 keV is expected in GRB 
spectra due to synchrotron self absorption, very high energy photons with 
£ > lO^^Ls2t~l_2T2\(ia/l keVr'eV still can escape from the GRB source 
(e.g. Li & Waxman (2007); see also Razzaque et al. (2004)). 
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assuming the > 100 MeV photons are produced in the same 
time and place as the MeV photons. Since thermal pressure 
acceleration in the initial fireball can not lead to much larger 
Lorentz factors, and in internal shock model much larger 
Lorentz factors would lead to synchrotron emission peaking 
below MeV band (see discussion in the first paragraph of §2), 
r2.5 « 1 is typically adopted (e.g Waxman 2003). The ex- 
act values of T would be determined by the detection of the 
high energy cutoffs in GRB spectra by, e.g. Fermi. We show 
below that the situation may be compUcated by the delayed 
large-radius emission from the outflow. 

3. LARGE-RADIUS EMISSION FROM RESIDUAL COLLISIONS 

After the initial strong collisions at small radii driven by 
the small timescale, ~ 1 ms, variation in the outflow, there are 
residual collisions continue to occur when the outflow is ex- 
panding to large radii. As the velocity fluctuation in the flow 
is being smoothed out by the on-going collisions, the delayed 
coUisions become weaker, and the postshock electron energy 
and magnetic field are smaller, which give rise to synchrotron 
emission at longer wavelengths. LW08 has well discussed the 
dynamics of the residual collisions and naturally explained the 
optical flashes from GRBs by this delayed residual emission. 
As mentioned in LW08, the energy density in the delayed col- 
lisions is dominated by the primary 7-ray emission, so that the 
residual emission is dominated by IC scattering of the prompt 
7-ray photons. In what foUows consider the IC emission. 

3.1. Dynamics 

Let us approximate the outflow by a sequence of 3> 1 in- 
dividual shells. LW08 has considered the simplified case: the 
shell masses are equal; the shells are initially separated by a 
fixed distance rfyar; the extent of the shells is much smaller 
than the radius of the outflow; and the Lorentz factors of in- 
dividual shells are drawn from a random distribution with an 
average T and initial variance cTp q < F^. The model may be 
more complicated, by adding more degrees of freedom, how- 
ever LW08 has shown that this simple case naturally accounts 
for the observed properties of GRB optical flashes. Here we 
wiU consider this simplified case. Moreover, we adopt the as- 
sumption that two shells merge into a bigger shell after a col- 
lision, i.e. the full inelastic collision where the internal energy 
generated is fully radiated. If the postshock electrons carry an 
energy close to equipartition and the electrons cool fast, the 
internal energy is always radiated significantly. In the case of 
a significant fraction of the internal energy being dissipated 
in each collision, the dynamical evolution of the outflow has 
been proven to be similar to the full inelastic case (LW08). 

In the simple case under consideration, the variance of 
the velocities of individual shells (in the outflow rest frame) 
evolves with the outflow radius R as cry oc Br^l^. The outflow 
energy that is associated with the fluctuation of shell veloci- 
ties (in the outflow rest frame) and hence may be dissipated 
decreases as 

£fluc cx oc (8) 

In general, it is naturaUy expected that there might be a wide 
range of variability timescale, ~ 1 ms - 10 s, in the flow prop- 
erties. Large timescale variabilities might lead to more energy 
dissipated at large radii. Thus, the slope should be flatter than 
-2/3. If a power-law description, Efluc oc is still available, 
one may have Q <q<2/2). We carry a monte carlo simulation 
to demonstrate this point in the appendix. 



3.2. High energy emission 

Based on the dynamical evolution, the emission from the 
residual collisions can be further predicted. Taking the com- 
mon assumptions that in internal shocks the postshock elec- 
trons and magnetic fields carry fixed fractions, and f^, re- 
spectively, of the postshock internal energy, the characteristic 
Lorentz factor of postshock electrons (in the outflow comov- 
ing frame) scales as 7,- oc e^o'l oc R~^l^, and the postshock mag- 
netic field scales as oc eBCyWe oc R~^l^ (the particle number 
density scales as Ue oc R~^). 

We demonstrate that the electrons in residual collisions lose 
most of their energy by IC cooling. If the prompt 7-rays last 
a duration T (observer frame), the plasma is overlapped 
with these 7-rays until the outflow expands to /? > 2F^cr ~ 
(T/tvarW-y- LW08 showed that when the synchrotron emis- 
sion lies in the optical band, the radius is /?opt — lO^R^. For 
typical observed values fvar ^ 10~^s and T ^ 10s, the op- 
tical radius is still relatively small, /?opt < 2r^cr. There- 
fore during the phase of late residual coUisions that we con- 
cern, the plasma is immersed in the radiation bath of the 
prompt 7-rays . Both the photon energy density {/- and 
the particle number density ng drop as oc R~^ hence the ra- 
tio y = Uj/(B'^/STr) oc oc increases with R. Because 
y ^ 1 in the 7-ray producing phase, we have y > 1 in resid- 
ual collision phase, so the radiation energy density dominates 
that of the magnetic field. Let us consider the properties of IC 
emission. 

3.2.1. Spectrum 

Consider first the energy band into which energy is radiated. 
At radius R, the prompt 7-ray photons with typical energy 
eh are up-scattered by electrons with characteristic Lorentz 
factor 7i to energy sic ^ X'yfeb \el(mp/mef{R/R^y^l'^eb. 
Here we assume 7,_o ~ Ceimp/me) as the electron Lorentz fac- 
tor emitting MeV 7-rays , and A accounts for the correction 
due to uncertain geometry effect. It will be shown in appendix 
that the correction factor A is order unity even in the case 
that the prompt MeV photons are strongly beamed in the rest 
frame of the outflow at R. The characteristic scattered photon 
energy is 

/ p \ -4/3 

The scattering might take place within slight Klein-Nishina 
regime, jic'ij/meC^ ~ a few > 1, at small radii R ^ R^, where 
the energy of scattered photons is instead eic = TjiMeC^ = 
eeT(R/R^)-^/^mpC^ ~ 3 x l0^e,(R/R^)-^/^GeV. 

Next consider the 77 absorption effect on the late residual 
emission. For GRB outflow with F > F^, the initial cutoff en- 
ergy for the primary emission is determined by the first term 
in the bracket of eq. (5), eo = ■ Eq- (3) impHes that the 
spectral cutoff energy scales as Scut RT'^L-\ so flie cutoff 
energy increases with R for fixed L. We have, for late residual 
collisions, ecut — soR/R-y. In the case of GRB outflow with 
lower Lorentz factor F < F^., the cutoff energy is initially a 

constant, £0 = £cut (^^■ 5), until the outflow expands to a ra- 
dius, 

R„ = 3x 10'^L52r2^5cm (10) 

(note/?m > R-y). At/? > the cutoff energy turns to increase 
with R, £cut — soR/Rm- In both cases of F > F^ and F < Fc 
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Fig. 1. — Schematic plot of the attenuated energy evolving with radius. 
There are two regimes. If the primary collisions that produce MeV emission 

occur in the regime of > R„ (i.e. F > Ft ), eg = £cut ^cut «: R later. 
If R-y < i?m (r < Fc), in the beginning the attenuated energy is a constant, 

£cut = £cut at R < Rm, and turns to be £cut <xRatR> Rm- 

the cutoff energy at /? > max[/?^, (i.e, Ecut > eo) follows 
the same expression, 

^ob / D \ 

.S.5„^(_)„.,r^GeV. ,U) 

The evolution of £cut versus R for fixed L is plotted in fig 1 . 

Comparing eic and £cut it can be found that typically eic > 
£cut at small radii, R < SOR^. In this case the bulk IC ra- 
diation is absorbed in the source, leading to electromagnetic 
cascades, and the photons escape until their energies decay 
to Ecut- Therefore the bulk high energy radiation is just re- 
emitted at Ecut- On the other hand, £ic < £cut at large radii 
R > 30Rj, where only the high energy ennission below £cut 
appears. The emission above Ecut is truncated, and undergoes 
electromagnetic cascades, but does not affect much the ap- 
parent spectrum since the photon spectrum decreases rapidly 
with energy. 

Finally consider the emission flux. It is straightforward to 
show that the cooUng time of the electrons is short compared 
to the dynamical time during the late residual collision phase, 
up to radii R ^ 10^ R^. We therefore assume that electrons 
radiate away all their energy. When eic > Scut at small radii, 
the total electron energy appears at Ecut — eoifR^<R<Rm or 
Ecut ocRif R> max[Rj,Rm]- The observed (time-integrated) 
IC spectrum at energy e > eq would be uF^, oc £'fluc|i/cu,=!^ oc 
^"^/'Ucu,=-^ oc (Here = e/h). 

When the outflow expands to large radii where Eic < Ecm, 
we need to consider the electrons accelerated to Lorentz fac- 
tors larger than the characteristic Lorentz factor 7,. Shock 
acceleration is expected to generate a power-law energy dis- 
tribution of electrons dne/d^e oc 7*'' at 7e > 7, with p ~ 2 
(Blandford & Eichler 1987). This flat-electron energy distri- 
bution, ^yjdrie/d^e oc 7^, generates equal amounts of IC en- 
ergy in logarithmic photon energy intervals, uF^, oc for 
ly > v\c- So vFyiy^ux) — i-'Fyiyic) fifiuc- The emission at 
low energy would be covered by earlier emission, while only 
the emission at high energy end, i.e. around the cutoff, shows 
up and interests us. The observed (time-integrated) spectrum 
would be similar to the eic > Ecut case, i.e., vF^, oc 



Thus we expect the observed (time-integrated) prompt 
spectrum, above the spectral cutoff energy in the prompt 7- 
ray emitting phase, eq (eq. 5), to be 

vF^ oc i/"^/', hu > £0 (simplified case). (12) 

This fluence spectrum is resulted from summing up all emis- 
sion components from different radii and times. A schematic 
plot of the prompt GRB spectrum is shown in fig. 2. Note, 
the spectral slope v'-^^^ here is derived from the simple case, 
which has been confirmed by recent numerical calculation by 
Aoi et al. (2009). If in general Efluc oc we would expect 
the slope to be uF^ oc 

Below Eo is the observed prompt MeV 7-ray spectrum, i.e., 
typically vF,, oc v below e^, and vF^ oc between Cb and eq. 
Note, the transition of the emission from primary to residual 
collisions at Eq is smooth if R^ > R,„, as shown by the dashed 
dot line. The transition for the case of R^ < R^ is discontinu- 
ous as shown by the thick solid line. The power law described 
by eq. (12) starts with a flux lower than the primary emission 
at £() by a factor of (Rm/R-y)^^^- 

Some comments should be made here. The spectral form 
described in eq. (12) holds only on average, especially for 
the high energy range. In individual GRB events the flux may 
differ significantly, because, for a small number of shells (and 
colUsions) at large radii, large variations in the late residual 
collisions should be expected. 

It should also be noticed that we have assumed the initial 
variance (Jr.o < r> whereas initial condition with trr.o > T 
may lead to more efficient 7-ray production (e.g. Beloborodov 
2000) around R^, in which case the ratio between fluxes of 
primary and residual emission at eq should be larger by a fac- 
tor of a few"*, leading to more abrupt transition between the 
primary and residual emission. 

One may worry about that the IC emission may be reduced 
as the seed photons are not isotropic in the shock frame of 
residual collisions. However, this kind of geometry effects 
do not play an important role even if the photons are com- 
pletely collinear (Wang et al. 2006). As usually assumed by 
many authors, suppose that the electrons accelerated in resid- 
ual collision shocks are isotropic in the rest frame of the out- 
flow, since the tangled magnetic fields in the shock might suf- 
ficiently isotropize electrons. Thus an electron is changing 
its angle 9' willi respect to the photon beam and cooling fast. 
The IC power of an electron averaged over its cooling time 
is not different from interacting with isotropic photons. As 
long as the jet effect is not important to prompt GRB emis- 
sion (^jet ^ 1 /r), we can furthermore regard the GRB explo- 
sion as isotropic, and hence the observer at different angles 
will observe the same IC emission due to spherical symmetry. 
Consider both cases of isotropic and anisotropic scatterings, 
the radiated energy can be assumed to be the same because 
it is determined by the total electron energy if electrons radi- 
ate all their energy rapidly. The outside observers would ob- 
serve the same fluence in both cases, otherwise one can sim- 
ply ask where the electron energy have gone, given the same 
total electron energy. Thus, the IC fluence is not reduced by 
this geometry effect if electrons are isotropic distributed in 
the rest frame and radiate all their energy within a dynamic 
time. Nevertheless, this effect changes the angular distribu- 
tion of IC emission. In the rest frame the IC power becomes 

It is not expected that the initial variance of Lorentz factors is far exceed- 
ing the mean, <Tr,o S> F. 
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P(q oc (1 -cos^^')^ although not much different from isotropic 
distribution. Correspondingly, in the lab frame the "image" 
of the IC emission is different from the isotropic case, i.e. 
the anglular dependence of the brightness is different, but the 
angular-integrated fluence is the same. 

3.2.2. Time delay 

At energy e < £o, the emission is mainly contributed by pri- 
mary collisions at small radius, and arrives at detectors simul- 
taneously with MeV emission. However, for higher energy 
£ > eo, the emission is produced at relatively large radii, and 
should have a time delay relative to the primary MeV emis- 
sion. Since the spectral cutoff energy, where the high energy 
photons emerge, increases as the outflow expands, Scut o<: R 
for Ecut > £o (eq. 11), the time delay increases with ob- 
served photon energy. For energy e = ecut(^) > £a, the typi- 
cal radius where photons emerge is ~ [e/eo] max[/?^, R„} = 
[e/eJ;[,{]/?-y, therefore the related time delay relative to MeV 
7-ray emission, Tdeiay — (R-R-y)/2T^c, is 

T^elay ^ — ^var = 2.1Ls2Tt,.^^^{l+z)h (for £* > sf). 

^cut 

(13) 

Thus the emission at e > eo is delayed later as e increases. 
A comment is made here that the target photons for scatter- 
ing may be beamed with respect to electrons, which changes 
the angular distribution of IC emission, i.e. the maximum 
IC power may come from a certain angle other than ^ = 0. 
This leads to an additional time delay. However, as implied 
by eq. (B5), the angle where the maximum power is emit- 
ted is smaller than 7r/2 in the comoving frame and hence 
d < 1/r. The produced time delay is smaller than Tdeiay. 
R6'^/2c < R/2T^c ~ Tdeiay Thus we neglect this additional 
time delay. 

Eq. (13) implies that the detection of time delay Tdeiay at e 
helps to determine the Lorentz factor. 

This determination by time delay should be consistent with 
that by detection of the spectral transition so between primary 
and residual emission, eq. (7). 

4. APPLICATION: FERMI-LAT GRB 080916C 

As the Fermi observational data showed up after the first 
version of this paper was posted on the archive,^ we now ap- 
ply the model to the observations. 

Abdo et al. (2009) reports the measurements of the bright 
GRB 080916c by Fermi GEM and LAT. The redshift of 
this burst is z = 4.35 (Greiner et al. 2009), which implies, 
with flucence (10 keV-10 GeV) w 2.4 x 10"'' erg cm 2, 
the largest reported isotropic 7-ray energy release, /iiso — 
9 X 10^"* erg. The observed GRB duration is T « 50 s, 
so the bolometric isotropic-equivalent luminosity is Lboi = 
£iso(l +z)/r ~ lO^'^ergs"'. As the observed peak energy is 
€i,/(l +z) ~ 1 MeV and the high energy slope is /3 w -2, 
the MeV luminosity, defined as the luminosity at < 2e/,, is 
L«Lboi/ln(10GeV/lMeV)= lO^^ergs'^ The LAT detected 
145 photons at > 100 MeV, within which 14 are beyond 1 

5 http://arxiv.org/abs/0810.2932 
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Fig. 2. — Schematic plot of the predicted i^/v spectrum of prompt high 
energy emission from a GRB. The thin solid line shows the observed MeV 
7-ray emission, a broken power law with a break energy at ej, ~ MeV, above 
which the spectrum goes as uF^ oc i^". The dot lines mark the break energy 
Cj, and the previously thought pair-production spectral cutoff, sq ~ GeV (eq. 
[5]). The residual collisions at large radii contribute beyond sq. Summing up 
all emission components (the dashed lines) from different radii and times lead 
to a spectral slope vFu oc v^i (q is the index of the random energy evolution, 
£fluc oc R^'', and < q < 2/3). In the simplified case, q = 2/3 (see the text). 
Here the thick solid line corresponds to the outflow satisfying < R„, (or 
£0 — TifieC^), while the thick dashed dot line corresponds to R-y > R„, (eo 
FmcC^). The uF,^ values in these two are different by a factor of {R,„/R-i9'^^ ■ 
The spectral slope holds up to TeV range, but not higher (see the discussion 
section). 

GeV, during the first 100 s after the trigger. The brightness 
may have enough statistics for spectral and temporal analysis 
of the high energy properties. 

There are several interesting properties in the high energy 
emission of this GRB. 

Time delay — The multi-band light curves unambiguously 
show that the bulk of the emission of the second light-curve 
peak is moving toward later times as the energy increases (see 
time bin b in Fig 1 and its inset panels in Abdo et al. 2009), 
and the time delay of lOO-MeV emission is about 1 s rela- 
tive to MeV emission, much larger than the MeV variability 
timescale, < 100 ms (Greiner etal. 2009).^ First of all, these 
are qualitatively consistent with our prediction that the higher 
energy photons can only arise when the plasma expands to 
larger size in later time where the 77 optical depth reduces to 
below unity, and that the size of high energy emission can be 
much larger than MeV emission. 

Let us consider the data quantitatively, and constrain the 
model parameters. The LAT >100 MeV detection consists 
of 145 photons which mainly come up in a single light-curve 
peak, therefore we have enough statistics for the time analysis 
of >100 MeV emission. It is obviously seen that the >100 
MeV light curve peak has a time delay r^i^^y = Tdeiay(l +z) ~ 
1 s relative to fliat of 250 keV-5 MeV (The script "ob" denotes 
quantities measured in the observer frame, with redshift effect 
taken into account). This implies £* = £o/(l +z) < 100 MeV. 
Substituting the observed values of L = lO^^ergs"', Tdeiay — 
1/5 .35 s and £ = 100 X 5.35 = 535 MeV, with redshift z = 4.35 

* Note, the time-delay issue here is different from what is called "delayed 
onset" by other authors. We concern indeed the delayed peaking time of high- 
energy emission, related to the delayed arrival of the bulk of high-energy 
emission. 
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taken, into eq. (13), we obtain the bulk Lorentz factor of GRB 
080916c outflow, 

/ J X 1/6 /^ob I \ "'/^ 

— (w4^) (^) ^ 

This result is similar to those determined in other GRBs, F » 
100-400, through observations of the rising of optical after- 
glows (Molinari et al. 2007; Kruhler et al. 2008, 2009; Greiner 
et al. 2008), and the thermal components in the prompt emis- 
sion (Pe'er et al. 2007). 

The determination of F can be double checked by the loca- 
tion of £o. Using the result of eq. (15), we obtain 

e^^""' ~ 10(?t/102ms)(T^,''i,y|iooMev/ls)-iMeV 

and 

e^^" ~ AQ{Ll\Q'\igs-'il\Ttx,y\mM.v/l s)-V«MeV. 

The observed MeV variability timescale is f < 100 ms based 
on the INTEGRAL observation (Greiner et al. 2009). Thus 
~ 40 MeV, consistent with requirement < 100 MeV. In 
addition, the broad Ught curve peak in the no-energy-selection 
band of LAT is consistent with, or a little delayed from, that 
of the GBM (260 keV-5 MeV) light curve, and is ahead of 
the >100 MeV peak. Thus ef should be located in the LAT 
energy window (no selection) and below 100 MeV, consistent 
with the result ~ 40 MeV. 

By our model the >1 GeV emission should be even 10 times 
longer delayed than the >1 00 MeV one, i.e., r^eiay 1 1 cev ~ 10 s. 
The much fewer photons above 1 GeV prevent us from ana- 
lyzing the temporal properties with high confidence. However 
the LAT >1 GeV Ught curve does agree with a longer delay 
by ~ 10 s. 

It should be noticed that other authors also constrain the 
bulk Lorentz factor of this GRB and obtain much larger lower 
limit (Greiner et al. 2009; Abdo et al. 2009). Essentially, the 
difference is due to different models considered; they con- 
sider the GeV emission produced in the same time and place 
as the MeV emission, whereas in our model the GeV emis- 
sion comes from delayed residual collisions at large radii, 
therefore our model looses the constraint on F. In addition, 
we consider that the cutoff energy should not locate below 
l^l'^TnieC^., which is ignored in Greiner et al. (2009) and Abdo 
et al. (2009). 

Lack of the first LAT light-curve peak — The low energy GBM 
light curve shows two peaks, however the LAT observations 
only show one peak related to the second GBM peak and 
there is a paucity of emission in the first ~ 4 s after the 
trigger. Note, some people call this as a "delayed onset" of 
high energy emission. In principle, one of the explanations 
could be that there is a spectral cutoff at ^ 10 MeV for the 
prompt emission from primary collisions in the first 4 s, and 
the residual-colUsion emission at >10 MeV is > 4 s delayed, 
longer than the second peak. If so, the properties of the ejecta 
emitting the first GBM peak are different from the later ejecta, 
which also suggests that there might be long-timescale, a-few- 
second i':^ ?var ~ 1 ms), variabiUties in the outflow of this 
GRB. 

Time-integrated spectrum — The joint GBM-LAT spectrum of 
GRB 080916c can be fit by Band function (Band et al. 1993), 



with peak energy around 1 MeV, a w -1 .0 and /3 « -2.2, ex- 
cept for the first 4 s (Abdo et al. 2009). Because the time 
intervals used to construct the spectra are much longer than 
the MeV variabiUty time, the resulted spectra are all time- 
integrated ones. Since the synchrotron self-Compton (SSC) 
model for GRBs, where the MeV peak is from the IC scat- 
tering by soft photon emitting electrons, predicts a bright 
GeV-TeV component due to the second order up-scattering 
(Piran et al. 2008), no evidence for high energy bump up to 
10 GeV in observations does not favor SSC but synchrotron 
model (Wang et al. 2009). In addition, the narrow vF^ spec- 
tral peak of GRB spectra favor more synchrotron emission 
mechanism over SSC, since SSC usually has a much broader 
spectral bump (Baring & Braby 2004). 

In the framework of synchrotron internal shocks, our resid- 
ual colUsion model predicts a slight spectral softening at high 
energies. However, due to the small detected GeV-photon 
number the Poisson scatter of low statistic still allows a slight 
softening at tens MeV to fit the data. Furthermore, if there are 
large timescale, ^ 1 ms, variabilities in the outflow so that the 
residual emission spectral slope q is larger than 2/3 (see §A), 
the high-energy spectrum is less steepened and is closer to a 
single power law. Finally, the slope of eq.(12) holds on aver- 
age, while the later residual collisions occur between smaller 
number of shells, thus there might be fluctuation from this 
average slope. 



Highest energy emission — The highest energy photon is de- 
tected wifli 13.2 GeV only 17 s after the GRB trigger. Wifli 
the redshift z = 4.35 this suggests GRB 080916c produces ra- 
diation up to 71 GeV in the source frame. Moreover, LAT 
detects 145 photons with > 100 MeV, within which 14 with 
> 1 GeV and especially only one with > 10 GeV, consis- 
tent with a power law spectrum with photon index /3 » -2 
up to ~ 10 GeV scale. There might be emission extending to 
higher energy, say, beyond tens of GeV, from GRB 080916c 
following the same slope, but the detection rate is less than 
one, i.e., no photon would be detected at this energy. Thus 
the observations actually suggest that the high energy spectral 
cutoff (or steep drop), if there exists, is more likely to be far 
above the energy of the only observed highest energy photon, 
^max = 13.2 X g GeV with ^ 3> 1. If the high energy emission 
beyond 13.2 GeV is produced by internal shocks, we argue 
here that it may not be produced by primary internal collisions 
that emit MeV 7-rays but produced in other regions, e.g., by 
residual collisions. 

As said above no high energy spectral component in GRB 
080916c does not favor SSC model but synchrotron model. 
Now calculate the maximum synchrotron photon energy. If 
B is the magnetic field strength in the internal shock, the Lar- 
mor time of an electron with Lorentz factor 7 is ?^ = ^rrieC / eB. 
The typical particle acceleration time can be scaled by Lar- 
mor time as = /f^ (e.g. Hillas 1984) where / is a correction 
factor accounting for the uncertainty of shock acceleration. It 
might be that / > a few (e.g. Lemoine & Revenu 2006). In the 
same time the electron suffers synchrotron cooling in a typical 
timescale f' = 3>mec/AaT^{B^ /^tt) (We neglect the IC cool- 
ing as the IC scattering usually occurs in deep Klein-Nishina 
regime for the most energetic electrons). The competition be- 
tween acceleration and coohng results in a maximum Lorentz 
factor of accelerated electrons, 7niax = (67re//c77-B)^/^. The 
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relevant synchrotron photon frequency is a constant/ 
^Lx = 0-'^lLxeB/2nm,c = 0.9e^ / fcfTm.c. 

The coefficient accounts for the fact that the synchrotron 
power per unit frequency of an electron peaks at the frequency 
a fraction 0.3 of the common characteristic frequency. The 
maximum synchrotron photon energy s'^^ = hiy'^^J^ / (I + z) 
is, therefore, 

e^ = 15r2.5rkl+zr'GeV. (16) 

Comparing the predicted maximum synchrotron energy 
with that implied by the observation, s'^^ = 13. 2g GeV, we 
have a lower bound, 

r=1.5xloVo.5,?o.5, (17) 

where the conservative values, / = 10*''''/o,5 and g = I0'^'^g{).5, 
have been taken. This bulk Lorentz factor is too large for 
fireball-shock model, because it faces several problems. First, 
the large F leads to (primary) internal shock radius larger than 
the deceleration radius of GRB outflow (e.g. Lazzati et al. 
1999). The deceleration radius is Ra « iEk/4nmpC^T^y^^. If 
Rj < Rd, an upper limit is obtained, 

r < 7 X lo3(£,,55/no)'/'?,t^!2> (18) 

where = 10^^£i:,55erg and n = l«ocm~^ are the outflow ki- 
netic energy and medium density, respectively. Second, the 
large T raises problem of low energy conversion efficiency 
due to slow cooling of accelerated electrons (e.g., Derishev et 
al. 2001). If the synchrotron cooling time of electrons with 
typical postshock Lorentz factor 7^ ~ mp/nie ~ 10^73, is re- 
quired to be smaller than the dynamical time of the outflow, 
t'cilm) <td^ Rj/Tc, we have 

r < 5 X 10^]/\Lboi/10^^Vgs-')'/^C/^2. (19) 

In this calculation we have assumed that the postshock mag- 
netic field is limited by observed emission, B^/Stt <U^ = 
iboi/47r/?^r^c. Third, in synchrotron internal shock mod- 
els the large T leads to large colhsion radius R^ « IT^cty^, 
where the magnetic field B is too small to give rise high en- 
ergy synchrotron photon energy. Using « TTij^^eB / nieC and 
the limit B^/Sn < Uj, the restriction to obtain synchrotron 
emission peaking at MeV range is 

r < 0.4 X 10373(Lboi/105Vgs-i)i/\-i/!2(ej,/lMeV)-V2. 

(20) 

Finally, the thermal pressure of the initial fireball is not ex- 
pected to accelerate the loaded baryons to very large Lorentz 
factor with most energy kept as the kinetic energy of baryons. 
The final Lorentz factor is limited to be 

r < 3 X 10^(Lo/10^''erg s'^^^'j^^ (21) 

where Lq is the rate at which the central source emits energy, 
and ro = 10^ro,7 cm is the source size (see, e.g Waxman 2003). 

' This upper bound for synchrotron energy is robust for any accelera- 
tion mechanisms involving electromagnetic processes, because the acceler- 
ation limit with / = 1 is robust not only to Fermi shock accelerations but 
also to any particle accelerations through electromagnetic processes. There- 
fore this boimd might be valid not only to internal shock models but also to 
electromagnetic-dominated models. 



The above upper limits to F imply that it may be impossible 
that F ^ 10^. This appear not to match the value suggested 
by the highest energy band observation, F = 1 .5 x 10Vo.5^o.5> 
unless / « 1 and ^ « 1 are satisfied at the same time: the 
shock acceleration must operate at the Bohm limit; further- 
more the observed highest energy photon happens to be at the 
maximum synchrotron energy. Thus the observations imply 
there is emission much higher than 13.2 GeV, which cannot be 
originated from synchrotron emission in the primary internal 
shocks. Actually, this high energy emission can be produced 
by IC emission in residual collision shocks, as discussed in 
present study. So the observation of highest energy emission 
supports the residual colhsion model. 

In conclusions: (1) the time delay of high energy emission 
and the spectral feature of highest energy emission in GRB 
080916c might have provide evidences for the residual colli- 
sion model; (2) its spectrum is not inconsistent with the resid- 
ual emission; (3) the time delay of > 100 MeV emission con- 
strains the bulk Lorentz factor to be F ~ 300, a typical value 
usually taken. It appears to be an applicable method to deter- 
mine F of GRB outflows by measuring the time delays of LAT 
light curves. If internal shocks also work in short GRBs, we 
expect similar delayed, prompt high energy emission in short 
GRBs. We also caution more careful spectral analysis to find 
the transition between primary and residual emission. 

5. DISCUSSION 

We have considered the high energy emission in the prompt 
GRB spectrum, which is dominated by the IC emission from 
the residual collisions. Instead of a exponential spectral 
cutoff, a steeper, compared to the prompt MeV emission, 
power-law slope i^F^, oc i^"* is expected beyond the previ- 
ously thought cutoff, £0 (eq. [5]). Here q is corresponding to 
the dynamical evolution of the random energy in the outflow, 
Efluc oc In the simplified case (see §3.1), which is consis- 
tent with optical flash observations (LW08), we take q = 2/3, 
while < (7 < 2/3 in general. The extended emission makes 
it complicated to detect the "cutoff energy" in the goal to con- 
strain the GRB emission region. 

Indeed, EGRET had detected prompt high energy photons 
past GeV in several brightest BATSE GRBs occurring in its 
field of view (e.g., GRB 930131, Sommer et al. (1994); GRB 
940217, Hurley et al. (1994)), which suggest that the other 
faint GRBs may produce prompt GeV photons as well (Din- 
gus 1995). There is also no sign of cutoff in the spectra (Din- 
gus 1995), which, if there is, should be far exceeding ~ 1 GeV. 
These EGRET results are consistent with the predicted ex- 
tension of prompt emission beyond GeV. However, the cut- 
off is not ruled out. Given the sensitive dependence of the 
cutoff on F (eq.5), a slight variation of F in individual GRB 
events may lead to much higher cutoff energy, » 1 GeV, ex- 
plaining the prompt GeV emission in EGRET-detected GRBs. 
Two properties may help to discriminate our residual-collision 
emission model from a very high energy cutoff model. The 
first is the steepening turnover in the spectrum. For typical 
Lorentz factors, F « 10^- 10^, the expected spectral turnover 
is ~ 100 MeV — 1 TeV, weU located in the windows of Fermi 
and AGILE. The second is the time delay of high energy emis- 
sion. One may expect systematic time delay of the high en- 
ergy photons in the residual emission model, while no delay 
is expected in the very high energy cutoff model. However, 
the task is not easy given that for a typical event with fluence 
^ 10"* erg cm~^ the observed GeV photon number is only a 
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few. In order to have enough statistics for the spectral and 
temporal analyses, very bright events are needed, or one may 
integrate many events to obtain an average burst. 

We have discussed that the recent Fermi detected bright 
GRB 080916c might have presented a good sample. Observa- 
tions do not show a simple spectral cutoff, but a spectral tail 
up to 70 GeV in GRB frame. More than one hundred of pho- 
tons detected above 100 MeV makes it obviously showing a 
time delay about 1 sec, which can be explained by the resid- 
ual IC emission and results in the determination of a typical 
Lorentz factor value, T ~ 300. The features of GRB 080916c 
support the residual-collision emission model as opposed to 
the very-high-energy-cutoff model. 

The high energy emission would not extend to very high 
energy. There are several effects that lead to a drop in TeV 
range. First, when the plasma expands to very large ra- 
dius, R > IT^cT - 10'^r^ 5(r/10 s) cm {T is the MeV 7- 
ray duration), there would be no overlapping between the 
plasma and the MeV 7-rays , and hence no scattering is ex- 
pected. At /? ~ IV'^cT the cutoff energy increases to £cut — 
l(r/10 s)f-i^_2 TeV (from eq.ll). A lack of > 1 TeV photons 
would be expected for GRBs with duration T < 10fvar.-2 s, 
although it should not be an exponential cutoff. Second, in 
interaction with ~ 1 MeV photons, the Klein-Nishina limit 
becomes important for electrons with 7,, > F, giving rise to 
IC photons up to eic — O.lFjjTeV, beyond which the spec- 
trum gradually turns below the low energy slope. Finally the 
cosmic infrared background would absorb the > 0. 1 TeV pho- 



tons that arrive from GRBs at redshift z ^ 0.5. Except for low 
redshift events, the observed prompt GRB emission may not 
extend far beyond TeV range. 

It should be commented here that in the framework of the 
synchrotron internal shock model (Waxman 2003) we do not 
expect a bright high energy component, say > 1 GeV, in the 
prompt emission, other than the synchrotron self-Compton 
model (e.g. Piran et al. 2008), therefore the residual high en- 
ergy emission will be dominant. A detection of high energy 
component in the prompt emission will be an evidence against 
the synchrotron model for GRBs, and vice verse. The recent 
Fermi observations of several GRBs do not support a high en- 
ergy component in GRB spectrum, since the > 100 MeV flu- 
ences are all less than those in MeV range (Abdo et al. 2009). 

Our residual collision model is not expected to produce 
much longer delayed high energy emission which is not ap- 
parently overlapped with the primary MeV 7-rays in times. 
There are more and more observations showing that GRBs 
produce delayed > 100 MeV emission even after the prompt 
7-rays end, lasting tens or even '--^ 10"* seconds (Hurley et al. 
1994; Gonzalez et al. 2003; Giuliani et al. 2008). This may 
require some long-lasting central engine activities or external 
productions (e.g. Wang et al. 2006, and references therein). 
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APPENDIX 



DYNAMICS WITH MULTI-TIMESCALE VARIABILITIES 



In order to show the effect of multi-timescale variabilities on the dynamics of the outflow, we carry simulations for both cases 
of single- and multi-timescale variabihties for comparison. We consider a series of individual material shells i = 1,2. ..,N, with 
total shell number A'' = 3000, released in a duration of 7" = 3s, so that the interval between two nearby shells is 1 ms. The shells 
have equal masses but different energies, with the bulk Lorentz factor of each shell following 

logr, = 2 + eiog9 + <^,-, (Al) 

where ^, is a random number between zero and unity and can be taken as the following forms, 

(Single); 

E,sin(^)logAi (Multi). ^^^^ 

For single-timescale case, we should take ^, = (case S), then the outflow has only one variabiUty timescale of 10"^s. The Lorentz 
factors randomly and uniformly distributed in logarithmical scale between 100 and 900. For multi-timescale effect we consider 
3 extra timescales besides the smallest timescale: = 10"^, 10"' , 1 s for k= 1,2,3, with relevant values being A^ = 2, 1.5, 1.2 
(case Ml) or A;t = 1.2, 1.5,2 (case M2) for k= 1,2,3. A^; decreases with in case Ml but increase in case M2, which means case 
M2 has larger fluctuations at larger timescales. Larger timescale fluctuations tend to produce strong collisions at larger radii so 
that larger fraction of energy is dissipated at larger radii. We further consider that in each two-shell colUsion, 1 /3 of the generated 
internal energy is emitted by radiation, because only the energy of shocked electrons is assumed to be emitted rapidly and the 
electron equipartition parameter is e^. = 1/3. The two shells are considered to separate again after collision and share equally the 
remained internal energy (in the center-of-momentum frame of the two shells). 

In Fig 3 we show the fraction of emission energy £'em(> ^)/£em,tot, that is emitted beyond radius R. We see that case S show a 
slope close to the analytical resolution oc for single-timescale case by LW08. However, the multi timescales lead to flatter 
slopes, implying more fraction of emission energy tends to be emitted at larger radii. Case M2 has even flatter slope than case 
Ml, since case M2 has relatively larger fraction of energy that is dissipated at larger radii. The steep drop at the end in both cases 
Ml and M2 means no more strong collisions later on. This is because there is a largest timescale of 1 s in our simulations. If 
there is still variabilities with timescale larger than 1 s then the slope will continue to even larger radii and show even later drop at 
the end. In summary, the simulations demonstrate that multi-timescale variabilities lead to a flatter slope q < 2/3, and the ^ = 2/3 
slope is only for single-timescale case. 
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Fig. 3. — The fraction of emission energy as function of radius in three simulations. The line marked with squares: case S, with only the 1-ms variability; 
circles. The other two are multi-timescale cases (Ml andM2), with three more scales ofPk = 10~^, 10~' and 1 s, for ^= 1,2, and 3, respectively. The line marked 
with circles: case Ml, with Aj. = 2, 1.5, and 1.2 for k = 1,2, and 3 (eq.[A2]); triangles: case M2, with = 1.2, 1.5, and 2 for = 1,2, and 3 (eq.[A2]). The 
straight line show the -2/3 slop to guide eyes. See the relevant text for more details. 



ANISOTROPIC IC EMISSION 

Let us discuss at which energy the IC emission is emitted in the observer frame, taking into account the fact that the seed 
photons are beamed. Consider the extremely anisotropic case, where the photons are coUinear in the comoving frame of the 
outflow. In this frame the electrons, as argued, is reasonably assumed to be isotropicaUy distributed. For simpUcity, we consider 
mono-energetic photons, since the photon number rapidly decreases with energy. Thus the IC power per unity solid angle in the 
comoving frame is angular dependent, 

gocd-^y, (Bl) 

where /x' = cos 6' with 6' the direction with respect to photon beam, and we have taken the velocity of the electron to be « 1 . 
Hereinafter prime denotes quantities in the comoving frame of the outflow, while non-prime denotes observer frame. Using the 
Lorentz transformation, /x' = (/i-/3r)/(l -/?rM). where /?r = (1 - 1 /r^)'/'^, we have 

l-M' = (l+/3r)-r^^, (B2) 
l-/3rM 

then the angular distribution of IC power in observer frame is 

dl7 r4(l-/3rA<)3df^' °^ (l-/3r/^)5' ' 
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This is not like the simple cone-like distribution of the isotropic-photon case. The maximum power per solid angle is emitted at 
angle with 

Mmax = ^^=^. (B4) 

The corresponding angle in the comoving frame is given by 

1 ' 2(l+/3r) ^4 

^-^™-=^:0^"5- ^^^^ 

The scattered photon energy in the comoving frame is given by e[Q « 7t?e'(l - fJ-'), with je the electron Lorentz factor and e' the 
photon energy both in the comoving frame. The photon energy (in observer frame) emitted at angel Umax, where the IC power is 
maximum, is then 

/ \ £ic(Mmax) 6 2-n /ri , a \ ^ 2 /•■d^:\ 

where in the last equality e = (1 + firWe' is taken for coUinear photons. This is the observer-frame photon energy around which 
the IC emission is mainly emitted. We see that the anisotropic correction factor A in eq. (9) is only order unity, implying that 
taking A ~ 1 is a good approximation even for highly beamed seed photons. 
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